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A careful study is made on the lattice of the phase diagram of QCD with two staggered 
flavors, to investigate the order of the chiral transition of Nf = 2 QCD. The specific 
heat and the susceptibility of the chiral condensate are determined for different spatial 
sizes of the system, and a finite size scaling analysis provides a determination of the 
(pseudo)critical indices. The result is a strong indication that the chiral transition is 
first order. 
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1. Introduction - Motivation 

QCD with Nf = 2 is a specially interesting system to understand confinement. A 
schematic view of the phase diagram is shown in Fig. ^ where for simplicity the 
two quark masses have been put equal to m and /z is the barion chemical potential.. 
At /X = the transition line between the region which is conventionally named 
"confined" and the "deconfined" one is defined by the maximum of a number of 
susceptibilities (Cy, Xjpipi ■ ■ ■) which happen to coincide within errors. As m — > 00 
the quenched case is recovered, and the transition is known to be first order, with 
(L), the Polyakov loop, as order parameter, and Z3 as symmetry. In principle Z3 is 
explicitly broken by the coupling to the quarks and (L) is not an order parameter: 
however the quenched description is valid empirically down to m ~ 2.5 — 3 GeV. At 
m ~ a chiral phase transition takes place, from the low temperature phase where 
chiral symmetry is spontaneously broken to a phase in which it is restored: there 
(il'tp) is a good order parameter. The mass term explicitly breaks chiral symmetry 
but (ipip) is expected to work as an order parameter in a neighborhood of m = . 
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Fig. 1. Schematic phase diagram of Nf = 2 QCD. 



At intermediate values of m neither (L) nor (t/'V'} expected to be good order 
parameters, even if their susceptibihties show a maximum along the transition line 
of Fig.n At m ~ another phase transition is expected to take place, with restora- 
tion of the Ua{^) symmetry, which is broken by the anomaly at low temperatures. 
An independent definition of confinement is needed to ask the question whether 
deconfinement transition coincides with chiral and C/yi(l) transition. An effective 
critical free energy density can be written for the chiral transitiorl^, assuming 
that the relevant critical excitation are the scalar and pseudoscalar particles. 

^ : 0,, = (9.(l+75)gj) ihJ = l-Nf) (1) 

Under chiral and ?7a(1) transformations of the group C/a(1) ® SU{Nf) SU{Nf), 
d> transforms as 



(2) 



so that by the usual symmetry arguments, and neglecting irrelevant terms 

- ^Tr [d.^^d^cf] - '^Tr{4>^} ~ y ffi [Tr{4>H}\ ' - ^g^Tr{4>^Y. (3) 

Inclusion of anomaly brings in an additional term L'^ — c \det(f) + det(f^^^ , which is 
SU{Nf) ® SU{Nf) invariant, but not Ua(1) invariant. 

One can inquire, by use of renormalization group plus e-expansion techniques, if 
infrared stable fixed points exist, which indicate the possible existence of a second 
or higher order transition. For Nf > 3 no such point exists, and the chiral transition 
is first order. At Nf ~ 2, in the absence of anomaly (c = 0) or if the r/ mass vanishes 
at Tc, no fixed point exists, the transition is first order, and also at m, /i ^ the 
transition is expected to be first order. In such a case no tricritical point exists in 
the plane (/x, T) of Fig.Qj^. If instead m^' ^ Q at Tc, and the Ua{^) transition occurs 
aX T > Tc, the symmetry group is 0(4) and the transition can be second order. In 
that case at m, /i 7^ there is no transition but only a crossover, and a tricritical 
point is expected in the (/u, T) plane of Fig. 2] 

The issue has fundamental implications for confinement. If the dcconfining tran- 
sition is order-disorder, so that an order parameter exists, a crossover is excluded. If 
instead a crossover exists, a state of free quarks can continously be deformed to the 
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confined phase, and can exist also there. The existing 

literaturJSEEll 

is not con- 
clusive on this point, even if a preference is given to the second order plus crossover 
scenario. The problem deserves more attention. We are going to present preliminary 
results based on a 5 TFlops x month lattice simulations on APEmille computer. 



2. Lattice Investigation 

The order of the chiral transition, and, more generally, the transition line of 
Fig. n can be studied by lattice simulations and standard finite size scaling tech- 
niques. We have used staggered fermions on Lt x lattices, with Lt = 4 and 
Ls = 12, 16,20,24,32. The input parameters of the simulations are (3 = 6/17^ and 
am, the quark mass in units of the inverse lattice spacing a^^. The temperature is 
given by 

T = ^- f4) 

with a{/3,m) the lattice spacing in physical units. The reduced temperature r = 
(1 - T/Tc) is then given by 

a[p, m) 

or, in a sufficiently small neighborhood of the critical point 
with 

In a 

(7) 



dm 



^ In g 



(/3c,0) 



df3 

In the quenched case fc = and t oc {(3c — P) ■ 

The correlation lenght of the order parameter ^ diverges at the critical point 

S^^T-" for T^0+ (8) 

with a critical index known as v. For the specific heat Cy and for the susceptibility 
X of the order parameter, the following scaling laws are expected 

Cv-Co = L':'-^c{TLl''',amqLyj^) (9) 
X-Xo^LZ/'^<^^{TL\/\am,Ly-) (10) 

Eas. l9llDl are obtained by renormalization group arguments, holding when the lattice 
spacing is much smaller than 1^, so that a/^ « 0. This is true for second order and 
weak first order transitions. The critical indices a, 7, yu depend on the transition 
(see Table ^ . The physics behind all that is that around Tc the free energy can 
be written in terms of the order parameter, and the dependence is dictated by 
symmetry: irrelevant (higher dimensional) terms are neglected. Eq|5|for the specific 
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Fig. 2. Specific heat (top) and Xm (bottom) peak value for Runl (left) and for Run2 (right), 
divided by the appropriate powers of Ls (^Eas. [TTll2l to give a constant. Horizontal stripes indicate 
the l(j confidence region for a fit with a constant value, /dof is also shown. 

heat is always valid. EqEH is only valid if the choice of the order parameter is 
correct. Therefore the specific heat is a reference to test the validity of any order 
parameter. 

In order to test the 0(4) option and second order, we have run at fixed value 
of the scaling variable amgif'' with — 2.49, and different spatial sizes Lg. The 
expectation is then, from Eqs. 191 101 that the peak values scale as 

{Cv - Co)^'""'= oc L^J" (11) 
(X - Xo)^^'^' « Lllr (12) 

Since the value of a for 0(4) is negative the height of the peak for Cv — Cq is 
expected to decrease at high Ls- Data show instead a rapid raise. The quantities 
{Cv - Co)P'=°V-^-s and (x - XaY^"^ I " . which should be constant if 0(4) were 
the symmetry, are displayed in Fig. [3 The /dof for a constant is very high: 0(4) 
symmetry is excluded. Since yh for 0(2) is equal within errors to that of 0(4), also 
0(2) symmetry can be tested and it turns out to be equally bad. Our action is not 
"improved": however we do not expect that an infrared property like Eas. llllT^ are 
affected by ultraviolet improvement. We can safely state that 0(4) and 0(2) are 
excluded and with them the crossover scenario. 



Table 1. Critical exponents. 
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Continuity arguments and Eas. lyillH require that, at small values of m 

Cy - Co ~ (am,)-"/(''^'')<i>c(TLy^ am,Lf ) (13) 
X-Xo- (am,)-^/(''^'•)$^(rLl/^ am,Lf ) (14) 
The positions of the peaks (3max scale then as 

Pc-Pma. + km-k'L;'^-' ^0 (15) 

their heights as 

{Cv - Co)P""''(am,)"/(''^'') = const (16) 

ix - Xof'' {amgy/'-'y^'^ = const (17) 

Eas. 1131 ini 1151 116l andll7lcan be tested with the data. Eq.^jis co mpat ible both 
with first order {ly = 1/3) and 0(4), in agreement with previous analyse d I I Eas. 1161 - 
Elagaiii exclude 0(4) and 0(2) and are consistent with a first order transition (see 
Fig El). 
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Fig. 3. Peak of the specific heat divided by the appropriate power of the mass f Eg. 1161 to give a 
constant, for different scahng hypotheses. 



3. Dual Superconductivity of the Vacuum 

A disorder parameter (fi) detecting dual superconductivity of the vacuum has been 
developed by our group during the last years and proved to be a good parameter for 
the quenched theory. The parameter can be defined equally well in the presence of 
dynamical quarks. {^) is strictly zero in the deconfine phase of Fig. ^ and it drops 
to zero at the transition line, as shown by the fact that the susceptibility 

P^IjMm) (18) 
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has a sharp negative peak coincident within errors with the peak of Cy . Around 

(/i)=i^><i>(^)(riy^am,L^'•) (19) 

continuity arguments require 

{^,)^{am,)-'^/^'^y-^^^,^{TLl''^) (20) 

and therefore 

p/Ll''^ :^f^^^{TLl'n (21) 

Independence on m is expected. Eq. [2] allows a determination of v. If this agrees 
with the value obtained from C^,, a legitimation of (/i) as order parameter follows. 
Eq-EJis compatible with first order (Fig.0J. 




Fig. 4. p/La scaling according to Eg. l21l for different scaling liypotheses. 



4. Conclusions 

Finite size scaling analysis of Nf = 2 QCD definetely excludes 0(4) -and 0(2)- 
symmetry at the chiral critical point, and favours a first order transition. The in- 
vestigation is continuing with larger lattices and improved actions, to better test 
the first order option. Dual superconductivity of the vacuum is confirmed as a good 
symmetry for the order parameter. 
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